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Abstract 

We investigate a (l+l)-dimensional nonlinear field theoretic model with 
the field potential V{cl)) = \(p\. It can be obtained as the universal small 
amplitude limit in a class of models with potentials which are symmetrically 
V-shaped at their minima, or as a continuum limit of certain mechanical 
system with infinite number of degrees of freedom. We point out that the 
model has an interesting scaling symmetry. One consequence of that is the 
existence of self-similar solutions of the corresponding field equations. We 
also find discontinuous solutions of shock wave type. 



PACS: 03.50.-z, ll.lO.Kk 
Preprint TPJU - 8/2005 



1 



1 Introduction 



Classical (l+l)-dimensional field theoretic models with nonlinear field equations 
appear in many problems in physics. Well-known examples are provided by sine- 
Gordon and 0^ models, which have plenty of applications: from Josephson junc- 
tions to macroscopic mechanical systems, see, e. g., fl]|2|. Very rich contents of 
these models have been explored in great detail. In the ultralocal limit, i. e., when 
gradient terms in the field energies are neglected, both models can be regarded 
as describing infinite sets of simple uncoupled one dimensional mechanical sys- 
tems: pendulums or anharmonic oscillators, respectively ^ In the present paper 
we consider a model which in that limit describes an infinite set of coupled balls 
bouncing vertically from a floor, but of course this is not the reason for investigat- 
ing it in detail. The actual motivation is that the model possesses very interesting 
scaling symmetry. Furthermore, it universally describes the small amplitude sec- 
tor in models with non smooth field potentials. In a sense, it is a counterpart of 
free field models which describe the small amplitude sectors in the case of smooth 
field potentials. 

Specifically, we consider the (l-i-l)-dimensional, classical field-theoretic model 
with the field potential 

^(0) = 101, (1) 

and the corresponding field equation 

^^ = -sign(</)(e,r)), (2) 

where is a real scalar field. The sign function has the values ±1 when 7^ 0, 
and if = 0. The potential V reaches its absolute minimum at </> = 0, where it 
is not differentiable - the left and right derivatives of V have different values: 

dV{(t)) 
d(f) 

Hence, the potential is symetrically V-shaped at the minimum. This feature dis- 
tinguishes our model from the majority of models of condensed matter physics 
and particle physics. Notice also that Eq.(2) can not be linearized even if is 
arbitrarily small. 

'Let us notice that the ultralocal limit yields systems with infinite number of degrees of free- 
dom. It should not be confused with the infinite wave length limit which leads to systems with one 
degree of freedom: a single pendulum or anharmonic oscillator. 



±1. 
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The non-smooth field potential (1) appears in an effective field theoretic de- 
scription (in a continuum approximation) of a discrete mechanical system com- 
posed of balls, which are connected by elastic strings and bounce vertically from 
the floor, see Fig. 1. Another system with a V-shaped potential - elastically cou- 
pled pendulums bouncing between two rods - is described in 01 . Moreover, V- 
shaped field potentials appear in theoretical investigations of pinning phenomenon 
of extended objects, see, e. g., H. These facts prove that eventhough the model 
might look rather exotic, it is not unrelated to the physical world. 

The model is quite intriguing. First, it is distinguished by the above mentioned 
universality. Second, apart from the obious l-i-l dimensional Poincare invariance, 
it possesses the exact scaling symmetry of the 'on shell' type, which resembles a 
little bit the approximate scaling symmetry of Navier-Stokes equations for a fluid 
in the turbulent regime. This symmetry is the most profound feature of the model. 
Because of its presence we may expect that there exist so called self-similar so- 
lutions of the field equation, see, e.g., |5| for the relevant mathematical concepts. 
Indeed, we have found such solutions. They have quite interesting structure both 
from the physical and mathematical viewpoints. Another very interesting feature 
of the field equation (2) is that it admits shock waves: finite discontinuous steps 
of the field which move with the velocities ±1 ('the velocities of light') and 
leave wakes behind them. One more reason for our interest in that model is that, 
to the best of our knowledge, effects of the nonlinearity of the sign(0) type in field 
equation have not been thoroughly investigated yet. 

The plan of our paper is as follows. In Section 2 two derivations of the model 
with potential (1) are given: (a) the continuum approximation to the system of 
bouncing balls, (b) the small amplitude approximation to a general symmetrically 
l^-shaped potential. Section 3 is devoted to the scaling invariance and to the self- 
similar solutions. The shock waves are discussed in Section 4. Section 5 contains 
a summary and certain remarks. 

2 The derivation of the model 

2.1 The continuum limit in the system of elastically coupled, 
vertically bouncing balls 

Balls of mass m can move along vertical poles which are fastened to the floor 
at the points = ai, i = 0, ±1, ±2, . . . , lying along a straight line. The balls 
bounce elastically from the floor. Each of them is connected with its nearest neigh- 
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bors by straightlinear segments of elastic strings. The system has one degree of 
freedom per ball. It is represented by the elevation h{xi, t) above the floor of the 
ball at the point Xi and the time t, see Fig. 1 . 




Figure 1: The system of balls connected by strings. The balls can move along 
the vertical poles (the continuous vertical lines) without friction. They elastically 
bounce from the floor which is depicted as the continuous horizontal line. Here 
e ~ /i, ^ ~ X and r ~ t are the dimensionless variables introduced in the text. 

Let us forget for the moment about the floor. Then each ball is subject to 
the force of gravity (represented by the gravitational acceleration g), and to the 
elastic forces from the two attached strings. It is an elementary exercise to write 
the Newton equations of motion for the balls. Next, we pass to the dimensionless 
variables: 



r 




= K^h{xi,t), 



where k is the elasticity constant for the strings, is a unit of length along the 
poles. Finally, we take the continuum limit: 



m 

0, — - = constans. 
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As the result, we obtain the equations of motion for the system of balls in the form 



provided that e(^, r) > 0. This last condition is due to the fact that we have not 
taken into account the presence of the floor. 

In case e = at certain point ,^0 and time tq, Eq. (2) is replaced at that point 
by the elastic bouncing condition: 

de de 

^(^0, i"o) ^^(^0, To) when e(^o, tq) = 0. (4) 

We can remove this cumbersome condition with the help of a trick used earlier for 
the system of pendulums, 01511, and called 'the unfolding' of the model. In the 
unfolded model instead of the field e(^, r) > we have a new field e(^, r) which 
can take arbitrary real values. By assumption, the evolution equation for e has the 
form (2), 

-^t^ - -^t^ = -^ign(i(^, r)). (5) 

The field potential corresponding to the r.h.s. of Eq. (5) has the perfect symmetri- 
cally V-shaped form, namely 

Z(e) = |e|. (6) 

Because the r.h.s. of Eq. (5) is finite, the velocity de/dr is a continuous function 
of time T. On the other hand, the two second order derivatives of e can not be both 
continuous when e = because the r.h.s. of Eq.(5) is discontinuous. 

The original field e is related to e by the following formula (the folding trans- 
formation) 

e(e,r) = |e(e,r)|. (7) 

It is clear that Eq. (5) and formula (7) imply Eq. (3) for e when e 7^ 0. The 
elastic bouncing condition (4) follows from formula (7) and from the continuity 

of de/dr: 

de/dr = sign{e)de/dr, 

hence the sign of de/dr changes when e (and e) passes through 0. 

The derivation of Eq.(5) presented above applies also to the phenomenon of 
pinning of an elastic string (e.g., a vortex) by a rectilinear impurity. In that case 
the gravitational force is replaced by the pinning force, and the system of balls 
and springs by the string. 
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2.2 The small amplitude sector for the symmetrically V-shaped 
potential 

Let us consider a general field equation of the form 



di 



de 



(8) 



where the field potential V{ip) is V-shaped at the minimum located at = 0, see 
Fig. 2. 




Figure 2: Generic V-shaped potential and the piecewise linear approximation. 

We are interested in the small amplitude oscillations of the field il) around the 
minimum of the field potential V . Because ip ~ 0, we may replace V by its 
piecewise linear approximate form 

v{ij) ^ v'{o+) ip Q{ip) - \v'{o-)\ ip e{-ip), 

where © denotes the step function. Then, Eq. (8) is replaced by 



-v'{Q+) e(^) + \v'{Q-)\ e(-v^). 



(9) 



Notice that this equation remains nonlinear for arbitrarily small V' - there is no 
linear regime for the small oscillations! After rescaling the field ^, 

^(e,r) = V'(0+)r/(e,r), 

we obtain the following equation 



-sign(?7) + 



\v'{o-)\ 

V'{0+) 



1 Qi-n). 



(10) 



Here we have used the formula Q{V'{0+)ri) = ©(r^) which is valid because 
V'{0+) > 0. Intheparticularcaseof symmetrically V-shaped potential I V(0—) I = 
V'{0+) and we obtain Eq. (2). 

Example of physical system for which Eq. (2) can be applied only in the 
small amplitude limit can be found in It is composed of elastically coupled 
pendulums which bounce from a rod. 

3 The scale invariance and the self-similar solutions 

Equation (2) has the scaling symmetry: if (f){^, r) is a solution of it, then 

0,(e,r)iA2 0(|,^), (11) 
where A > is an arbitrary positive number, obeys Eq. (2) too ^ . The field energy 

scales as follows: 

E[<P^]=X'E[<P]. 

When A — s> the solutions eA(^, r) are in general characterized by high fre- 
quencies, short wavelengths and small energies. This is in contrast to, e.g., the 
massless 0'^ model which also has a scaling symmetry, but in that case the scaling 
points to low frequencies and large wavelengths when the energy E[(j)x\ — > 0. 
Let us also remark that this reminds the phenomenon of turbulence - even more 
so when we recall that the Navier-Stokes equations in the high Reynolds number 
regime have a scale invariance similar to the one described above, see, e. g., the 
Introduction section in Q. 

The immediate consequence of the scale invariance is the lack of a characteris- 
tic frequency or energy scale. For example, one can find infinite periodic running 
wave solutions to Eq. (2) which have the dispersion relation of the form 

2 1 2 2 

uj — k = /i , 

where //^ can take any value from the interval (0, oo), 0. 
^Also Eq. (10) has such symmetry. 
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Let us stress that the symmetry transformation (1 1) is defined only in the space 
of solutions of the field equation (2). The action functional corresponding to that 
equation, 

S[<P]= jdrdi [1(9.0)2 _ 1(5^0)2 

which can be calculated also for fields which do not obey Eq. (2), is not invariant 
with respect to the transformations (11). Therefore, the symmetry (11) is not of 
the Noether type. In such a case one often says that the symmetry is of the 'on 
shell' kind, see, e. g., [8|. 

As always when there is a symmetry, one can find solutions of the pertinent 
field equation which are invariant under the symmetry transformation. In order 
to find the self-similar solutions of the field equation (2), that is the ones which 
are invariant under the scale transformations (11), we adopt the following scale 
invariant Ansatz for the solution: 

In consequence, Eq. (2) is reduced to the ordinary differential equation for the 
function S{y) 

{l-y^)S" + 2yS' -2S = -%\gn{S). (12) 
This equation has the following polynomial solutions: 

1 cy 1 

When 5>0: S{y) = --^{y^ + I) + -y + 

1 n' 1 

When 5<0: S{y) = -^'{y' + I) - -y - 

where a, /?, a', (3' are arbitrary constants. It has also the trivial solution 

So{y) = 0. 

The polynomial solutions are valid on appropriate finite intervals of the y-axis 
which are determined by the conditions 5 > or S* < 0, respectively. 

Let us try to put together the polynomial solutions in order to cover as large 
as possible interval of the y axis. The resulting solution has the form depicted in 
Fig. 3, where Sk-, /c = 1, 2, . . . , denote the polynomial solutions. We assume that 
5(0) = because ?/ = corresponds to the spatial infinity, ^ = ±oo, and we 
would like to minimize the growth of 0(,^, r) in that region. 
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Figure 3: The piecewise polynomial solution, are the matching points for the 
consecutive polynomials Sk- 

Let us write the polynomial solutions in the form 

Equation (12) implies that the following matching conditions at the points y — 
have to be satisfied: 



'S'fc(afc) — — 5'fc+i(afc), Sl.{ak) — S'^_^_i{ak), 



(13) 



provided that ak 7^ ±1. The points y = ±1 are exceptional because then the 
coefficient in front of the second derivative term in Eq. (12) vanishes. This has 
the consequence that the first derivative dS/dy of the solution does not have to be 
continuous at the points y = ±1. Below we shall make use of this possibility. 
The matching conditions (13) give the recurrence relations: 



Aai 



1-4 



fc+i 



1-4 



2afc - (1 + (>'k)(ik-i 
1 + a| - 2akak-i 



where ai — ai > 0, Pi = 1, qq — 0. We also assume that ai < 1 in order to 
ensure that 02 > oi. Rather unexpectedly, we have found exact solutions of these 
relations - it has turned out that a glance at explicit forms of several initial ak, f3k 
and ak is sufficient in order to guess the solutions. They have the following form 



^(l+ai)q'-' (1 - q^'-') , (3, = ^(l-a^) {l + q'-') (l + q-') , (14) 

(15) 



4 
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where 

1 — fli 



Of course, this solution can be checked by substituting it into the recurrence rela- 
tions. 

Because < g < 1, formula (15) implies that < 1, and that — > 1 when 
k ^ oo. Therefore, the piecewise polynomial solution constructed above covers 
only the interval [0, 1). For this reason we introduce a special notation for it, 
namely Spp{y). Simple calculations show that Spp vanishes when y ^ 1. On the 
other hand, the first derivative Spp{y) does not vanish in that limit, but it remains 
finite. 

This solution can be extended to the full range of y. First, we extend it to 
the interval (—1, 1) by taking as the solution in the interval (—1, 0] the function 
—Spp{—y), which smoothly matches the function Spp{y) at the point y = 0. Next, 
we take the trivial solution So{y) = for y < — 1 and y > 1. This is possible 
because Spp{y) and —Spp{—y) vanish in the limit y — > 1. The first derivatives 
of the functions Spp{y) and So{y) at the points y = ±1 do not have to be equal 
because Eq. (12) allows for a finite jump of the first derivative S'{y) when = 1. 

The parameter oi remains free except for the restriction < ai < 1, so we 
actually obtain a family of self-similar solutions. A snapshot of the solution at cer- 
tain r > is depicted in Fig. 4. It is composed of infinitely many quadratic poly- 
nomials in ^ taken on domains which become smaller and smaller when ^ — * r+ 
or ^ ^ — r— . For ^ ±oo |0| approaches the function aiT\^\/2 — The 
zeros of lie at Cfc(T) = T/ak- Hence, = ^/a-k > 1 - the zeros move with 
'superluminal' velocities. The function r) exactly vanishes for |{| < r, i. e. 
inside the Tight cone'. This follows from the fact that the restriction — 1 < y < 1 
is equivalent to |^| > r. 

The solution presented in Fig. 4 is symmetric with respect to ^ — > — It 
turns out that the left- and right-hand halves of this solution taken separately are 
solutions too: (/)(^,t) = e{^-T)C^S{^) and(/)(^,T) = e{C+T)C^S{^) also obey 
Eq. (2). 

The full self-similar solution described above can not be physically realized 
in a physical system. The reason is that it has infinite energy - this is the general 
feature of nontrivial self-similar solutions. Nevertheless, finite pieces of it can 
be observed if the dynamics of the model is local, as in our case. Namely, if a 
physical finite energy initial configuration differs from our solution only at the 
very large values of |^| (|^| r), this difference will be seen in the region of 
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101 




Figure 4: The self-similar solution. The two arrows indicate that the waves em- 
anate from the points which move along the ^-axis: ^ = ±r. We plot |{/)| because 
this is the physical variable in the case of system of balls - see the relation (7). 

finite |^| (|^| ~ t) only after certain finite time interval. During that time interval 
the piece of the self-similar solution describes the evolution of the finite |^| part 
of the system very well. Information about uses of self-similar solutions can be 
found, e. g., in Q. 



4 The symmetric shock waves 

Another interesting feature of the system governed by Eq. (2) is that it can support 
shock waves. Such waves are not invariant under the scale transformation (11). 
The Ansatz 

0(e,r) = 0(-z)iy(z), z = -^{e-r^). (16) 
reduces Eq. (2) to the following ordinary differential equation 

zW" + W' = ?,\gn{W), (17) 
where ' denotes the derivative d/dz. We have used the following formulas: 

e\-z) = ~6{z), z6{z) = 0, z6'{z) = -6{z), 

sign{e{-z)W{z)) = e{-z)sign{W{z)). 
Equation (17) has the following partial solutions: 

when W > : W+{z) = z + zi + zi \n\ — \ + di, 



11 



when 1^ < : W-{z) — —z — Z2 — Z2 In | 



Z2 



d2: 



where di,d2,zi, Z2 are constants. These solutions are defined on finite intervals of 
the z-axis, see Fig. 5. 

W+ 



II 



Figure 5: The partial solutions W±. The shock wave is obtained by smooth match- 
ing of infinitely many copies of the parts /, //. 

Putting the solutions W± together we obtain the solution which is defined for 
all 2; < 0. Because of the step function in the Ansatz we do not need to know W{z) 
for 2; > 0. In order to analyse the matching conditions for the partial solutions we 
introduce the following notation: 

W^^{z) =^ + ao, Wk{z) = (-l)'^ f^ + afe + Zfeln^ ) , 

where ak > 0. The matching conditions at the points —ak, /c = 0, 1, . . . , have the 
form 

(a) W^,(-a,) = = iyfe_i(-a,), (b) W^-a^) ^Wl_^{-ak), 
see Fig. 6. 

The matching conditions (b) for the derivatives yield the following recurrence 
relation 

Zk+i = 2ajt+i - Zk, where Zq = 2ao. (18) 

Here Oq > is a free parameter. The matching conditions (a) for the functions 
Wk also give a recurrence relation, but it is much more complicated. In order to 
write it in a transparent form, instead of we equivalently use Xk defined by the 
following formula 

ajfc+i = ao xi X2 • . . . • Xfe+i, Xjt > 0, A; = 0,1,.... (19) 
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W{z) 




Figure 6: The scheme for matching the parts /, /J of the partial solutions W±. 
Then, conditions (a) are equivalent to the following infinite set of equations for 

Xi, X2, ■ ■ ■'■ 

Xi = 1 + 2 Inxi, 
X2 = 1 + 2(1- lnx2, 

; (20) 

Xfc = 1 + 2 fl - ^ + i . . . ± i ) Inxfc, 

These equations can be rewritten in terms of Lambert W function [9|, but we 
will not give here the related formulas. Numerical solutions of the first twenty 
equations (20) are presented in Fig. 7, while Fig. 8 shows the corresponding 
values of a^. 

Xfi 
4 r 

3.5 • 

3 
2.5 

2 - . 
1.5 ' . ,^ 

1 

0.5 

2 4 6 8 10 12 14 16 18 20 " 
Figure 7: The numerical solutions x„ of Eqs. (20). 
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a„ 
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80 
60 
40 
20 



2 4 6 8 10 12 14 16 18 20 

Figure 8: The numerical results for the coefficients a„. 

The numerical calculations convincingly indicate that a;„ — 1 as n increases. 
Unfortunately, this is not sufficient to determine whether a„, related to by for- 
mula (19), increase to infinity or converge to a finite value. We have numerically 
investigated the sums J2n=i upto N — 250. The results suggest that the sum 
becomes divergent when N ^ oo. Therefore, we conjecture that a„ grow to infin- 
ity when n — oo. In this case the function W{z) is defined in the whole half-axis 
z <0, and formula (16) gives the solution for all values of z. 

The snapshots of the shock wave at three times are shown in Figs. 9, 10, 11. 







\v\ 




-To 


To 



Figure 9: The symmetric shock wave at an initial time tq such that < 
To < 2^ya^. The arrows indicate the directions in which the wave fronts move. 

(/.(e,To) = Ofor|e|>ro. 
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Figure 10: The symmetric shock wave at a time ti > tq such that 2.,Ja{ > t\> 
2^. Now Tl) = for \i\ > Tl. 

101 




-T2 




T2 



Figure 11: The symmetric shock wave at a time T2 > Ti. At later times more 
zeros of cj) appear. 

The features of the shock wave shown on Figs. 9-11 can be obtained as 
follows. The condition z <0 implies that r) vanishes outside the 'light cone', 
i. e., when > r^. The Tight cone' corresponds to 2; = 0, where = oq. The 
constant > is a free parameter. It gives the height of jump of the value of the 
field on the 'light cone'. Thus, Q{—z)W{z) describes a shock wave, symmetric 
with respect to ^ — > — ^ and restricted to the 'light-cone'. The velocities of the 
steps (shock fronts) are equal to ±1. 

At the given time r > the physically visible range of z is finite - it is given 
by the inequalities — < 2; < 0. Therefore, we can see only a finite number 
of zeros of 0, but that number grows indefinitely with time. The positions of the 
zeros of the field (f) are determined from the equations z = — 0^. Hence, they are 
given by the formulas 

Their absolute velocities are greater than 1. Thus, the zeros move behind the front 
and faster than it, but they never catch up with the front. The first zero of (p appears 
at the point ^ = at the time = 2y^. Notice that it is by factor ^/2 larger 
than the time a/Soq needed by a freely falling ball in order to hit the floor if its 
initial velocity is equal to zero and the initial elevation above the floor to Oq. The 
reason for this difference is obvious: the ball located at { = 0, which indeed has 



15 



the elevation oq and the vanishing velocity in the limit r ^ 0, does not fall freely 
because it is coupled to its neighbours. 

The discontinuous shock waves have infinite energy because the gradient en- 
ergy density at the shock fronts is infinite. Nevertheless, similarly as the self- 
similar solutions, they can be very helpful if used with due care. 

5 Summary and remarks 

The nonlinear term — sign(0) in the evolution equation (2) rarely is considered 
in field theory. Nevertheless, it is quite interesting, in particular because it has 
the scaling symmetry (11), and in consequence allows for the self-similar waves. 
Furthermore, we have shown in an earlier paper, see |3 1, that in other models with 
nonlinearity of that kind there exist static kinks which have strictly finite extension 
(no exponential tails) - for that reason they are called compactons. 

Systems with non-smooth potentials are quite common in classical mechan- 
ics. Especially interesting are so called impacting systems, for example, a har- 
monically oscillating particle which can hit a wall when it moves sufficiently far 
from its equilibrium position. Recently, such systems have been investigated in 
connection with chaotic behaviour and so called grazing bifurcation, see, e. g., 
ifTOlfTTl . In the present paper we have considered the field-theoretic counterpart 
of such systems. 

There are many interesting possibilities for future investigations of the models 
with V-shaped field potentials. Let us point two of them: 

1 . It is not known whether the models of that kind can be integrable. 

2. What are properties of quantum version of the (l-i-l)-dimensional model with 
the potential V{(f)) = |0|? We expect that the scaling symmetry (1 1) is broken in 
the quantum theory. What is the resulting mass scale? 
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